Entangling the optical frequency comb: simultaneous generation of multiple 2x2 and 
2x3 continuous-variable cluster states in a single optical parametric oscillator 
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We report on our research effort to generate large-scale multipartite optical-mode entanglement 
using as few physical resources as possible. We have previously shown that cluster- and GHZ-type 
iV-partite continuous-variable entanglement can be obtained in an optical resonator that contains 
a suitably designed second-order nonlinear optical medium, pumped by at most 0(N 2 ) fields. In 
this paper, we show that the frequency comb of such a resonator can be entangled into an arbitrary 
number of independent 2x2 and 2x3 continuous- variable cluster states by a single optical parametric 
oscillator pumped by just a few optical modes. 



INTRODUCTION 

Continuous-variable (CV) quantum information [l[ 
presents interesting distinctive features, first and fore- 
most its well developed experimental interface by use of 
quantum optics, where the field amplitude and phase 
quadratures are the exact mathematical analogs of po- 
sition and momentum of a quantum harmonic oscilla- 
tor. In particular, the Bloch-Messiah reduction Q pro- 
vides an elegant method for generating any iV-partitc 
entangled CV quantum state having a Gaussian Wigner 
function by the interference of N field modes emitted by 
N independent single-mode squeezers. Examples have 
been given for CV Greenberger-Horne-Zeilinger (GHZ) 
states Q and, more recently, for CV cluster states 0, HI- 

Cluster states, originally described using qubits @, 0| 
but later extended to qudits Q and continuous vari- 
ables Q , are the core resource for one-way quantum com- 
puting, a paradigm of quantum computation in which the 
unitary evolution of a quantum algorithm is implemented 
via measurements on a highly entangled state of a par- 
ticular type (a cluster state) , in a similar way to the use 
of quantum teleportation as a quantum computing prim- 
itive 11 1 . In addition, cluster states of any dimen- 
sion possess remarkable entanglement resilience to mea- 
surement: when one measures a subsystem of a cluster 
state in an appropriate basis, the rest of the system stays 
entangled in a cluster state. Cluster states are locally 
equivalent to Bell states for two entangled systems and 
to GHZ states for three entangled systems but become 
locally nonequivalcnt to those for four or more entangled 
systems. Recently, Grover's algorithm was implemented 
on fo ur q ubits using an optical one-way quantum com- 
puter [lj|. 

A one-way CV quantum computer based on Gaussian 
cluster states can perform universal quantum computing 
provided that at least one non-Gaussian detection oper- 



ation is used [^, [l3| . Such a detection operation could be 
a photon-number resolving stage, towards which signifi- 
cant progress has recently been made 1J| , with the other 
stages being straightforward high-fidelity homodyne de- 
tection. Therefore, striving toward the large-scale quan- 
tum engineering of Gaussian states seems a particularly 
worthwhile endeavor. 

In recent work ID, 3], we have shown that the afore- 
mentioned N independent squeezers and 0{N 2 ) beam- 
splitters can be replaced with a single .ZV-mode squeezer, 
based on a nonlinear photonic crystal (a quasi-phasc- 
matched periodically poled |20| ferroelectric) 

pumped by 0(N 2 ) field modes, and no output interfer- 
ometer. We outlined detailed experimental proposals for 
the generation of CV tri- and quadripartite GHZ 15|, 21 1 
and cluster [l6[ states. In these works, the entangled 
systems are different eigenmodes of the optical paramet- 
ric oscillator (OPO) cavity, labeled by their resonant 
frequencies and polarizations. Such a physical system 
exhibits high classical coherence, as spectacularly illus- 
trated by the use of optical frequency combs consisting 
of ~ 10 6 locked oscillating modes of a femtosecond laser 
as a scale- universal time/frequency standard 22,23(. R- e_ 
cently, ultrastable phase-locked OPOs were used to ob- 
serve macroscopic Hong-Ou-Mandel interference [24[ and 
bipartite CV entanglement [SB]]. The question we ask 
here is, can the optical frequency comb be turned into a 
quantum computer register? 

The theoretical answer is yes [TBI, EH- However, while 
there is no fundamental limit to the generalization of CV 
multipartite entanglement to the whole frequency comb, 
the experimental implementation does face one hurdle: 
the consequences of unwanted nonlinear interactions of 
modes within the cluster with others outside the cluster. 
These consequences arc difficult to assess precisely be- 
cause of the absence of good multipartite entanglement 
measures, but can be easily seen to be detrimental. We 
therefore adopt the conservative approach of designing 
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the iV-mode squeezer perfectly — i.e., with only the de- 
sired couplings between the modes inside the entangled 
set. The use of the selection rules provided by addi- 
tional quantum numbers, such as polarization, is a log- 
ical avenue to solve this problem. However, the two- 
dimensional polarization basis does not offer enough de- 
grees of freedom to scale experimental implementations 
beyond 4-mode entanglement. Nevertheless, we show in 
this paper that the simultaneous generation of multiple, 
larger CV cluster states is possible in the frequency comb 
of a single OPO, using only a few pump fields and simple 
nonlinear media. 



MULTIPARTITE ENTANGLEMENT IN THE 
OPTICAL FREQUENCY COMB: 
EXPERIMENTAL OUTLOOK 

In quantum optics, the set of resonant modes of an op- 
tical cavity, analog to the harmonics of a vibrating string, 
can be seen as a set of quantum harmonic oscillators 
whose frequencies are equally spaced, in the absence of 
dispersion, by the free spectral range (a "harmonic quan- 
tum harp" ) . If a nonlinear medium is inserted in the cav- 
ity and pumped by a monochromatic mode at frequency 
w pnmp , pairs of modes (m,ri) of respective frequencies 
that u> m +Lu n — Wpump (temporal phase-matching condi- 
tion, or photon energy conservation) will become coupled 
(Fig. Q}. This coupling manifests itself by the creation or 
annihilation of a photon pair with one photon in mode 
m and one photon in mode n. It is well-known that such 
an OPO yields bipartite entanglement below 2f| 27 1 and 
above [II, 0, the OPO threshold. The principle of 
the method lies in considering whether an appropriate 
network of such bipartite couplings may lead to multi- 
partite entanglement. As already mentioned, solutions 
always exist for GHZ and cluster states, although they 
are far from being unique (TBI . fl(| . In the following, we 
give a detailed description of the method. 



The aforementioned Hamiltonian is of the form 

H = ifay^G M («L«l - a m a n ), (1) 



where k is an overall coupling strength, a m is the photon 
annihilation operator for cavity mode m, and G = (G mn ) 
is a square symmetric matrix describing the Hamiltonian 
coupling network: we make the hypothesis of equal cou- 
pling constants for all coupled modes, so that the ele- 
ments of G are either or 1. (We will relax this later to 
allow —1 as well.) The matrix G is thus the adjacency 
matrix of a graph representing 7i, which we refer to as 
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FIG. 1: Physical system and corresponding TL-graph for a sin- 
gle pump mode. The flow of time is from bottom to top of the 
figure. The horizontal direction is the optical frequency axis 
for the signal OPO modes, with the pump modes denoted by 
half their frequency. The Ti-graph is drawn below and corre- 
sponds to mere pairwise entanglement. 

the 7i-graph as a reminder of this correspondence. An Ti- 
graph has vertices denoting the field modes (in a vacuum 
state) and edges representing the non-zero terms of H. 
The 7i-graph allows us to easily visualize the couplings 
associated with the implementation of the Hamiltonian 
in Eq. ([T]) . We sometimes refer to states created by inter- 
actions of the form of Eq. |T]) as 7i-graph states, or even 
simply 7i-graphs if this will not lead to confusion with 
the actual graph itself. Fig. [1] displays the example of 
a monochromatically pumped (type-I) OPO, which can 
only generate entangled pairs. The 7i-graph adjacency 
matrix in this case, denoted G\, has a constant main 
skew diagonal: 



G l = 



/ 1 \ 

10 

10 

y i o o o j 



(2) 



In our recent work [la ] , G was constrained to be sym- 
metric, real, and full rank. In this paper, we give more 
emphasis to the experimental viewpoint and add an- 
other constraint. Since the optical modes are labeled 
by their frequencies, the temporal phase-matching rela- 



H (Hamiltonian)-graph states: physical description tion (a->„ 



lp ) together with the assumption 



of constant interaction strength yield an adjacency ma- 
trix G with constant skew-diagonals, that is, with ones at 
m + n = (const.) for each given pump mode. Any matrix 
with constant skew-diagonals is called a Hankel matrix. 
For the sake of experimental simplicity, we therefore re- 
strict the discussion in this paper to G"s of Hankel form, 
even though more general ones can be considered (and 
may also be feasible) flo) . 

Throughout the rest of the paper, we adopt a more 
compact notation for Hankel matrices, in which we only 
list the entries along the top and down the right side of 
the matrix, with the top-right entry set off with slashes: 



Gi = [0,0,0/1/0,0,0]. 



(3) 
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This shorthand provides a way to read off the experi- 
mental implementation immediately in the case of a cav- 
ity with an evenly spaced frequency comb of modes (and 
no polarization dependence): each nonzero entry in this 
"vector" corresponds to a pump beam. The numerical 
value of the entry is proportional to the amplitude of 
the pump (negative would indicate a 7r-phase shift), and 
pumps corresponding to adjacent entries are apart in fre- 
quency by twice the free spectral range of the OPO cav- 
ity. 



Square-cluster OPO 

Unlike qubit grap h states (which include square-lattice 
cluster states) [30|, the entanglement properties of TC- 
graphs have not yet been studied in detail, except for 
our recent existence proof of an anal ytic correspondence 
between 7i-graph and cluster states [lfjj. It is, however, 
remarkable that very simple 7i-graphs on the frequency 
comb may yield fairly elaborate cluster states. Consider 
the physical system of Fig. [21 where two pump modes 
interact in the nonlinear crystal to yield four interacting 
modes. Note that it is crucial that the red connections 
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FIG. 2: Physical system and corresponding Ti-graph for a bi- 
raodal pump. Here, the interactions connect four modes each. 
As is shown in the text, these form a square cluster state. 

be limited to only the one drawn in the figure. For ex- 
ample, if there were a fifth mode in the cavity, modes 2 
and 5 (not drawn) must not be coupled by the pump 
mode drawn in red. This can be achieved by engineering 
the nonlinear optical interactions in quasi-phase-matched 
materials I la. Il9j . such as periodically poled KTiOPC>4 
(PPKTP) [3l|. Referring to Figf2] (and remembering that 
pumps are drawn at half their true frequency), we can 
write down the compact form of the adjacency matrix by 
inspection: 



G 2 = [0,0,0/1/0,1,0] 



(4) 



Under these conditions, it can then be shown that the 
4-mode set actually generates a CV cluster state with a 
square graph. More precisely, the defining relation [jjj 
for CV cluster states is satisfied, 



where p and x are column vectors of quadratures for 
each field mode, A is the adjacency matrix for a (possi- 
bly weighted) graph, and the arrow denotes the limit of 
large squeezin g,. N ote that the source of entanglement is 
a single OPO [16| . Moreover, the present experimental 
proposal is simpler than that of Ref. . 

A general proof of the above correspondence is based 
on the aforementioned analytical link which we will see 
in the next section. We now give a particular-case proof 
of our claim that the 7i-graph created by G 2 corresponds 
to a CV cluster state. Solving the Heisenberg equations 
for mode evolution under Eqs. ([T]) and (dJ), we obtain the 
following squeezed joint operators: 



(-S-X! + 5-X 2 - X 3 
(S + Xi - S + X 2 + X3 

(6-I\ - 

(5+Pi 



X 4 )e 
X A )e 



—rS. 



+ 



6+P 2 -Px + P 



4 e 



— r5- 
—r5+ 
r5— 



(0) 
(7) 
(8) 
(9) 



where r = nt and 5± = (v5± l)/2. Next, we compare 
these to the cluster-state equation, p — Ax — > 0, using 
the adjacency matrix A given by 



.4 



/ 0-2 -y/5 \ 

y/5 2 

-2 V5 

V -y/5 2 / 



(10) 



which corresponds to a weighted square graph with a 
certain choice of weightings. This yields 



P l + 2X3 + V5X 4 
P 2 - VdX 3 - 2X A 
P 3 + 2X t - V5X 2 
Pi + VlX 1 - 2X 2 



(11) 
(12) 

(13) 

(14) 



p -Ax — > 0, 



(5) 



which can be easily seen to be equivalent to Eqs. (JBHHl) , 
in the limit r — * 0, if we rotate modes 3 and 4 by ir/2 
(Ih>F, — X). This completes the proof. 



Multiple-square-cluster OPO 

Principle and first experimental implementation 

As seen in Fig. [IJ a single pump mode can generate an 
arbitrary number of entangled pairs in the same OPO, at 
least within the phase-matching bandwidth, which can 
be made orders of magnitude larger than the OPO's free 
spectral range. This somewhat trivial fact can be gener- 
alized to the square cluster of Figj2j with the same num- 
ber of pump modes (two), an arbitrary number of square 
clusters can actually be generated, as depicted in Fig. [3] 
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FIG. 3: Physical system and corresponding Ti-graph for a bi- 
raodal pump. Here, the interactions define three (or more) 
sets of four modes each. As is shown in the text, each of 
these form a square cluster state. 

for three copies. This corresponds to a Hankel matrix of 
the same pattern as G 2 but for AN modes rather than 
4. Once again, we can immediately read off the Hankcl- 
shorthand form of G3 by inspection of the diagram: 

G 3 = [On/1/05,1, 5 ] , (15) 

where 0„ is shorthand for a string of n zeros. From an ex- 
perimental perspective, the few pump modes and single 
type of nonlinear interaction greatly simplify the setup. 
One still has to prevent the red pump from introducing 
spurious mode couplings (i.e., couplings of the indicated 
modes to other modes outside the set), for example by 
designing sharp quasiphasematching cutoffs in the non- 
linear crystal or by filtering mode amplitudes or engineer- 
ing dispersion in the OPO cavity. Such a constraint isn't 
needed for the green pump because any spurious cou- 
plings generated would be strictly outside of the desired 
mode set. 

Second experimental implementation 

Note that this experimental generation of multiple 
square clusters is far from unique and can be even more 
efficiently and conveniently implemented by use of the 
polarization degree of freedom, as depicted in Fig. [4] 
This yields several advantages: first, one uses but a single 
pump mode; second, one has the advantage of connecting 
entangled modes in closed sets, without any additional 
engineering of the nonlinear interaction bandwidths. The 
interactions yield exactly the same equations as Eqs. (|llf - 

m. 

The experimental implementation is relatively 
straightforward: one can use one ZZZ crystal placed in 
sequence with a YZY crystal (first letter denotes pump 
polarization and other two the signal polarizations), with 
one crystal rotated at 90 degrees from the other. In this 
case, a vertically polarized pump generates the necessary 
VHV , VVH, and VVV interactions. Moreover, if the 
crystals have the same length, the fact they are rotated 
by 90 degrees from each other ensures that the frequency 
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FIG. 4: Physical system and corresponding H-graph for a sin- 
gle pump mode. The letters H and V denote the optical polar- 
izations of each mode, which are frequency degenerate. Here, 
the interactions define three (or more) sets of four modes, 
each of which form a square cluster state. 

combs corresponding to each polarization will have the 
same free spectral range, which in turns ensures common 
resonance conditions at a given cavity length. 



Third experimental implementation 

A variant of the previous implementation allows for 
generation of multiple copies of the exact proposal of 
Ref. (l6| . which yields a more balanced square cluster 
state. This is achieved by adding an additional interac- 
tion (Fig. dashed lines) of opposite effect to the three 
former ones (i.e. upconverting if the three are downcon- 
verting, or vice versa). This can be ensured by seeding 
the OPO below threshold with coherent-state signals, ap- 
propriately phase-shifted with respect to the pump fields. 
The experimental implementation consists, for example, 
in using 2 PPKTP crystals: one that simultaneously 
quasiphasematches the Y ZY and ZZZ interactions [31| 
and one only the ZZZ interaction, rotated at 90 degrees 
from the other. In this case, vertically and horizontally 
polarized pumps generate interactions VHV and VVH , 
VVV, and the opposite HHH. 

The resulting lower weighting [l6| can be seen, as be- 
fore, by solving the quantum evolution for the multimodc 



squeezing, to be compared to Eqs.([6][9]). 

(X 1 +X 2 -V2X 4 ) ) e~ rV * (16) 

(x 1 -X 2 + V2X 3 ) e- r ^ (17) 

(Pi + P% + V2 P 4 ) e- r ^ (18) 

(p - P 2 - V2 P 3 ) e"^, (19) 



and deducing, as before, the following cluster state rela- 
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FIG. 5: Physical system and corresponding H-graph for two 
frequency degenerate pump modes. Here, the interactions de- 
fine three (or more) se ts o f four modes, each of which form 
a square cluster state \lh\l . Note that the opposite interac- 
tion (dashed line) is not necessary but gives a more balanced 
cluster state. 



tions, to be compared with Eqs. (fTTl - [T4[) : 
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These implementations are currently underway at the 
University of Virginia. 

These methods therefore yield an arbitrary number of 
disconnected 4- mode square cluster states from just a sin- 
gle OPO, limited only by the number of modes within the 
phase-matching bandwidth. The whole frequency comb 
can thus be partitioned into entangled subsets by us- 
ing only one or two pump modes. Our eventual goal, 
of course, is an arbitrarily large square-lattice grid (or 
other universal graph), and many disconnected squares 
isn't the same at all. Still, the fact that only one or 
two pumps are needed to generate so many copies of this 
(rather modest) cluster state leaves us hopeful that this 
is a worthwhile line of research to pursue. In the next 
section, we give the theoretical basis for this result and 
show that this process can also be used for 8-mode cubic 
cluster states. These can each be converted, if desired, 
into 2x3 square-lattice cluster states at the expense of 
homodyne detection on two modes per cube. 



MULTIPARTITE ENTANGLEMENT IN THE 
OPTICAL FREQUENCY COMB: THEORETICAL 
OUTLOOK 

Mathematical relation between H- and cluster 
graphs 

The relation between a bipartite CV cluster state, rep- 
resented by the adjacency matrix 



A 



A( 
Al 



(24) 



and the aforementioned 7i-graph matrix G is given 
by 



G 



[B-_A CAl] [BA + AoC] 
[A^BAo - C] 



CAl 



A%B] 



(25) 



where B and C are arbitrary symmetric positive-definite 
matrices. This relation allows for abundant (in fact, ex- 
cessive) freedom in the choice of G for a given A. For 
concreteness, let's take B = C = I, where / is the 



identity matrix. 

AnA~ 



0^6 



Now assume that A is unitary, i.e. 
I. Then Eq. ([25]) yields 



l A 
1 A% 



A . 



(26) 



Restricting to unitary A's therefore can be used to elim- 
inate the distinction between CV graphs and 7i-graphs, 
simplifying the problem greatly (at the expense of some 
generality). Keeping in mind the possible limitations of 
such a restriction, we focus only on unitary A's. Using 
= to indicate equality up to a renumbering of modes, 
we then inquire whether we can always find a G = A, 
where G is in Hankel form, for obvious experimental con- 
venience. 



Simultaneously generating multiple copies of a CV 
cluster state 



For reasons that will become clear shortly, let's be- 
gin with a bipartite graph having a unitary, Hankel Aq. 
Since Aq is Hankel, it is automatically symmetric. Thus, 
Eq. ([24]) can be expressed as A = F 2 <£> A , where 



F n = [0 n _i/l/0 n _ 



(27) 



is the n x n skew-identity. Since Aq is unitary, by the 
argument above, we can choose for our 7i-graph, 



G = A <g) F 2 = F 2 (E) A = A 



(28) 



since exchanging the order of a tensor product is equiva- 
lent to a renumbering of the modes. It's straightforward 
to show that since Aq is Hankel, then Aq (Ei F n (for any n) 
is also Hankel. Considering that F 2 m (g> Aq is zero except 
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for 2TV blocks of Aq on the main skew-diagonal, it can be 
seen that it is the adjacency matrix for TV distinct copies 
of the graph corresponding to A. Also, F 2N ® Aq is uni- 
tary if A is unitary. Since F 2N ® Aq = A g) F 2N , the 
latter being Hankel if Aq is Hankel, then for any unitary, 
Hankel Aq, we know that 



G N = A Q ® F 2N 



(29) 



is Hankel and will create TV copies of the CV graph A as 
defined in Eq. (f24]) . 



Simultaneous generation of 2 x 2 and 2x3 cluster 
states 



To apply these results, let's first consider the case of 
multiple square clusters in a single OPO. A unitary Aq 
for a single square cluster is [r| 

A ° " 71 (? = 7i M/1/11 ' (30) 

which is also a Hankel matrix. Using the results above, 
we have that 

G N = Aq ® F 2N = [o, -1, 0/1/0, 1,0], (31) 

where = Q2N-1, is a Hankel matrix for the 7i-graph 
that will generate TV disconnected squares from a single 
OPO using just three pumps. 

Note that the negative matrix elements in Gn are re- 
quired for unitarity but are absent in Eq. (|15[) in the ex- 
perimental section, which yields a CV cluster state with 
the same graph but different weights. As we previously 
indicated, requiring that A be unitary is not necessary; 
we do it here only for the convenience of the calcula- 
tions. It proves useful, though, because we can extend 
this method to generate an arbitrary number of discon- 
nected 8-mode cubic cluster states, each of which can be 
reduced by measurement to a 2 x 3 square-lattice cluster 
state, if desired. We show this presently. 

The unitary form of Aq for a cubic cluster is 



Aq 



-1,-1,1/0/1,1,-1]. 



(32) 



The corresponding graph is depicted in Fig. \§\ Note that 
a cubic cluster state reduces to a 2 x 3 one under the 
measurement of two neighboring vertices in the position 
(amplitude quadrature) basis. Using the results above, a 
generating 7i-graph matrix in Hankel form for TV discon- 
nected cubes is 



G N = Aq® F : 



2N 




O 



O 





FIG. 6: Cubic cluster graph. Dashed lines indicate negative- 
weight edges. Measurement of any two neighboring vertices 
(open circles) in the position (amplitude quadrature) basis dis- 
entangles those vertices from the rest of the graph, thereby 
'ing a 2 x 3 cluster. 



-^[0, -1,0, -1,0, 1,0/0/0, 1,0, 1,0, -1,0] , (33) 



where, again, = 02at-i. We stress that this genera- 
tion of an arbitrary number of disconnected cubes can 
be achieved using only six pumps and a single OPO. The 
maximum number of cubes is limited only by the number 
of modes within the phase-matching bandwidth. 



CONCLUSION 

In this paper, we have expanded our previous results 
[l6| on the generation of CV cluster states with a min- 
imum amount of physical resources, in particular by 
adding the constraint that the 7i-graph adjacency ma- 
trix be of Hankel form and unitary. Although justifiably 
limiting the generality of the treatment, this constraint 
makes it easier to connect to concrete and simple ex- 
perimental cases, for which we show that large sets of 
independent 2x2 and 2x3 square-grid cluster states 
can be obtained in a single OPO. Multiple single-cube 
cluster states can also be produced in this fashion. 

A remarkable feature of this sort of Hankcl-unitary du- 
plication procedure is that the number of pumps is inde- 
pendent of the number of squares or cubes being gener- 
ated, as is the case for multiple entangled pairs; we need 
only three pumps at most (and as few as one, as shown in 
the experimental sections) to generate multiple squares 
and six pumps to generate multiple cubes. This is sig- 
nificantly better than our earlier estimate of the 0(TV 2 ) 
pump fields needed to create a cluster with TV modes. 
This estimate was based on the conservative assumption 
of one pump per graph edge and TV(TV— l)/2 ~ TV 2 edges 
in a complete graph. We observe that square-lattice clus- 
ter states are significantly sparser, with O(N) edges. In 
addition, single pumps can couple many mode pairs when 
Hankel matrices are used. Work is currently underway to 
experimentally generate these CV cluster states and to 
explore what other types of states can be created under 
these conditions. 
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